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ABSTRACT 



Context. Over the last decade a series of results have lent support to the hypothesis of the existence of a long thin bar in the Milky 
Way with a half-length of 4.5 kpc and a position angle of around 45°. This is apparently a very different structure from the triaxial 
bulge of the Galaxy. 

Aims. In this paper, we analyse the stellar distribution in the inner 4 kpc of the Galaxy to see if there is clear evidence for two triaxial 
or barlike structures, or whether there is only one. 

Methods. By using the red-clump population as a tracer of the structure of the inner Galaxy we determine the apparent morphology 
of the inner Galaxy. Star counts from 2MASS are used to provide additional support for this analysis. 

Results. We show that there are two very different large-scale triaxial structures coexisting in the inner Galaxy: a long thin stellar 
bar constrained to the Galactic plane (|Z>| < 2°) with a position angle of 43?0±1°8, and a distinct triaxial bulge that extends to at least 
\b\ < 7.5° with a position angle of 12°6±3°2. The scale height of the bar source distribution is around 100 pc, whereas for the bulge 
the value of this parameter is five times larger. 

Key words. Galaxy: general — Galaxy: stellar content — Galaxy: structure — Infrared: stars 



1. Introduction 

There is now a substantial consensus for there being a stellar bar 
in the inner Galaxy. This was first suggested by de Vaucouleurs 
(1964) to explain the nonaxisymmetry in radio maps. The first 
evidence for a barlike distribution in the stars was derived from 
the asymmetries in the infrared (IR) surface brightness maps 
(e.g. Blitz & Spergel 1991; Dwek et al. 1995) and in source 
counts (Weinberg 1992; Hammersley et al. 1994; Stanek et al. 
1994), which both show systematically more stars at positive 
galactic longitudes I < 30 close to the Galactic plane (GP). The 
exact morphology of the inner Galaxy, however, is still contro- 
versial. While some authors refer to the bar as a fatter structure, 
around 2.5 kpc in length with a position angle of 15-30 degrees 
with respect to the Sun-Galactic Centre direction (Dwek et al. 
1995; Nikolaev & Weinberg 1997; Stanek et al. 1997; Binney 
et al. 1997; Freudenreich 1998; Lopez-Corredoira et al. 1999; 
Bissantz & Gerhard 2002; Babusiaux & Gilmore 2005), other 
researchers suggest that there is a long bar with a half length 
of 4 kpc and a position angle of around 45 degrees. It is note- 
worthy that those authors supporting the 23° bar all examine the 
region at |/| < 12°, whereas those supporting the long bar with 
the larger angle are trying to explain counts for 10° < I < 30° 
and-10°>Z>-30°. 

Hammersley et al. (2000) showed that there is a major over- 
density of K2-3III (red-clump) stars on the GP at I = 27° at a 
distance of about 6 kpc from the Sun. This over-density could 
also be detected at smaller galactic longitudes, but more red- 
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dened and at an increased distance from the observer. These stars 
could not be seen a few degrees above or below the GP at the 
same galactic longitudes, nor at greater longitudes. This, when 
combined with the findings of Hammersley et al. (1994), was in- 
terpreted as strong evidence for the existence of a long in-plane 
bar with a half length ~4 kpc and position angle of 43°, clearly 
distinguishable from the triaxial bulge, but running into it near 
1=12°. This over-density of stars was also analysed by compar- 
ing near-infrared (NIR) counts with predictions of the Besancon 
Galactic model (Robin et al. 2003) with similar conclusions 
(Picaud et al. 2003). This result has been further supported by 
observations in the mid-infrared with GLIMPSE data (Benjamin 
et al. 2005) who did a similar analysis to that of Hammersley et 
al. (2000) but at a wavelength range in which the effects due to 
extinction are even lower than in the NIR, hence reducing the 
uncertainty in the results. 



There is also a noticeable difference in the luminosity func- 
tion found in the two regions. The majority of the work on the 
short, fat bar (subtending an angle of ~23° with the Sun-Galactic 
Centre line), has been restricted to within a few degrees of the 
plane in regions that are routinely used by other authors to study 
the bulge, and what is seen is an older population. The on-plane 
region between I = 20° and 27°, however, has been shown 
to contain a large number of extremely luminous young stars 
(Hammersley et al. 1994; Gazon et al. 1997; Lopez-Corredoira 
et al. 1997), and a characteristic feature of long bars is a major 
star formation region where they interact with the disc. 
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These results suggest that the structure in the inner Galaxy 
(|/| < 12°) could be interpreted as a triaxial bulgeQ rather than 
a bar, and certainly not a long straight bar. This triaxial bulge 
has axial ratios of 1:0.5:0.4 (Lopez-Corredoira et al. 2005). 
Furthermore Lopez-Corredoira et al. (1999) showed that a tri- 
axial bulge + disc model cannot reproduce the observed counts 
in the Galactic plane for 10° < / < 30°. This is also sug- 
gested by NIR photometry of red-clump stars (Nishiyama et 
al. 2005). It should be noted that the possibility of there being 
even smaller non-axisymmetric structures (e.g. a secondary bar) 
in The Milky Way cannot be discounted (Unavane & Gilmore 
1998; Alard 2001; Nishiyama et al. 2005); thus, more observa- 
tions are needed for a decisive conclusion to be reached. 

Many galaxies have triaxial bulges contained within a pri- 
mary stellar bar, and even a third, smaller-scale bar inside the 
bulge (Fiedli et al. 1996). More than half of all disc galaxies 
have some kind of bar component (Buta, Crocker & Elmegreen 
1996) and as many as a third have a strong primary bar (Freeman 
1996). Recent NIR observations by Beaton et al. (2005) and 
iV-body simulations by Athanassoula & Beaton (2006) provide 
strong evidence in favour of the existence in our close neigh- 
bour M31 of a boxy bulge contained within a primary bar of 
semilength of ~4 kpc. Such a high degree of morphological 
similarity between two neighbouring galaxies is striking (even 
though the angular separation between the major axes of the pri- 
mary bar and triaxial bulge is different for the two galaxies). 

Bars are generally constrained to the plane whereas bulges 
are shorter and more extended vertically. Therefore, using data 
covering a wide range of longitudes and latitudes is the only way 
to distinguish between the two possible morphologies for the in- 
ner Milky Way: a short inner bar or a triaxial bulge + long bar. 
In this paper we have used NIR data from the TCS-CAIN sur- 
vey to analyse the spatial distribution of red-clump sources in 
the inner Galaxy. This population of stars has been widely used 
as a standard candle in deriving distances to the inner Galaxy 
(Hammersley et al 2000; Stanek et al. 1994; Lopez-Corredoira 
et al. 2004; Babusiaux & Gilmore 2005; Nishiyama et al. 2005, 
2006b). By observing how these stars are distributed along dif- 
ferent lines of sight towards the inner Galaxy we are able to ex- 
plore the 3D morphology of our Galaxy. 

2. TCS-CAIN catalogue 

TCS-CAIN catalogue is a NIR survey of the Milky Way that 
has recently been completed at the Instituto de Astroffsica de 
Canarias. This survey consists of more than 500 fields distributed 
along or near the Galactic plane (\b\ < 10°) observed in J, H and 
K s filters with a photometric accuracy of ~0. 1 mag in the three 
bands and a positional accuracy of ~0.2", based on the 2MASS 
catalogue as astrometric reference. Each field covers approxi- 
mately 0.07 deg 2 on the sky, centred on the stated galactic co- 
ordinates for that field. The nominal limiting magnitudes of the 
survey (determined using growth curves of the differential star 
counts) are 17, 16.5 and 15.2 in J, H and K s respectively; hence, 
TCS-CAIN is at least one magnitude deeper than 2MASS, and 
even more in the inner Galaxy, where confusion is important. A 
complete description of the catalogue and data reduction, and de- 
tails of its contents can be found in Cabrera-Lavers et al. (2006). 



1 In other papers (Lopez-Corredoira et al. 2000, 2005) we have pre- 
sented evidence that the Galactic bulge has a boxy morphology. The 
(/, b) range examined in the present paper, however, is too restricted to 
make a proper contribution to that debate; hence, for simplicity we refer 
here to a triaxial geometry. 



In this paper we have used all the available near plane fields 
in the inner Galaxy < 30°, |£>|<7.5°), which gives a total num- 
ber of 205 usable fields, covering an area of approximately 14.2 
deg 2 of the sky (nearly 35% of the whole catalogue area). 

3. Red-clump stars as a distance indicator 

Red-clump giants have long been proposed as standard candles. 
They have a very narrow luminosity function and constitute a 
compact and well-defined clump in an HR diagram, particularly 
in the infrared. Furthermore, as they are relatively luminous, they 
can be identified even to large distances from the Sun. Lopez- 
Corredoira et al. (2002) developed a method to obtain the star 
density and interstellar extinction along a line of sight by extract- 
ing the red-clump population from the NIR colour-magnitude 
diagrams (CMDs). This method has proved very powerful for 
analysing the structure of the thin (Lopez-Corredoira et al. 2002, 
2004) and thick (Cabrera-Lavers et al. 2005) discs of the Milky 
Way, as well as in the study of the distribution of interstellar ex- 
tinction (Drimmel et al. 2003; Picaud et al. 2003; Duran & van 
Kerkwijk 2006). 

The absolute magnitude (Mk) and intrinsic colour, (J - 
K s )q, of the red-clump giants are well established (Alves 
2000; Grocholski & Sarajedini 2002; Salaris & Girardi 2002; 
Pietrzynski et al. 2003). Here, we consider an absolute mag- 
nitude for the red clump population of -1.62 + 0.03 mag and 
and an intrinsic colour of (J - K^o = 0.7 mag. These values 
are consistent with the results derived by Alves (2000) from 
the Hipparcos red-clump and also with the results obtained by 
Grocholski & Sarajedini (2002) for open clusters. The intrinsic 
colour of the red-clump stars predicted by the Padova isochrones 
in the 2MASS system (Bonatto et al. 2004) for a 10 Gyr popu- 
lation of solar metallicity is (J - K) = 0.68 mag. As 2MASS 
photometry is nearly equivalent to the TCS-CAIN photometric 
system (Cabrera-Lavers et al. 2006) the value of 0.7 mag is suf- 
ficiently close to represent the red-clump population. 

The absolute magnitude of the red-clump has a small de- 
pendence on metallicity (Salaris & Girardi 2002). There is no 
large metallicity gradient in the Galactic disc (Ibata & Gilmore 
1995a,b; Sarajedini et al. 1995); however, there is a slight dif- 
ference between the solar neighborhood metallicity and those of 
the bulge giants, which peaks at [Fe/H] = -0.25 dex (Ramirez et 
al. 2000; Schultheis et al. 2003; Molla et al. 2000). According 
to Salaris & Girardi (2002), a small correction of 0.03 mag 
should be applied to the absolute magnitude of the red-clump 
stars for these ranges of metallicity. However, this is well within 
the uncertainty considered for the absolute magnitude used here. 
The red-clump absolute magnitude in / is more sensitive to 
metallicity and age than the K filter, hence the intrinsic colour 
(J - K)o also depends on both the metallicity and age (Salaris 
& Girardi 2002; Grocholski & Sarajedini 2002; Pietrzynski et 
al. 2003). The theoretical isochrones of Salaris & Girardi (2002) 
and Bonatto et al. (2004) suggest that the de-reddened colour 
for the range of metallicities expected for a bulge population is 
(J - K)q = 0.63 + 0.02. Hence the metallicity dependences lead 
to a systematic uncertainty of 0.03 mag in the absolute magni- 
tude of the red-clump stars and 0.05 mag in the intrinsic colour, 
(J - K) . 

A recent study using NIR spectra of the inner Galaxy giants 
between 20° < I < 30° has shown that these stars have a well 
defined metallicity (typically [Fe/H] between and -0.2 dex) 
with a very small dispersion (Cabrera-Lavers et al., in prepara- 
tion). By using the population of infrared carbon stars from the 
2MASS survey, Cole & Weinberg (2002) conclude that the bar 
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formed more recently than 3 Gyr ago and must be younger than 6 
Gyr. For this age and metallicity, the corrections given by Salaris 
& Girardi (2002) are even lower than for the older bulge. Hence, 
we have assumed that uncertainties in absolute magnitude and 
intrinsic colour for this population of red-clump stars are of the 
same order as in the case of the bulge stars. 

In this paper the TCS-CAIN survey data is analysed us- 
ing the method developed by Stanek et al. (1997) for their 
data obtained by the Optical Gravitational Lensing Experiment 
(OGLE; Udalski et al. 1993, 1994). This method has also recently 
been used by Nishiyama et al. (2005, 2006b) and Babusiaux & 
Gilmore (2005) to analyse the red-clump population of the inner 
Galaxy. 

To avoid extinction effects a reddening-independent magni- 
tude (K e ) is used: 



1=7° b=2° 



K e — K s 



A,- A 



-( J - K s ), 



(1) 



adopting A K JEj- Kx - 0.68 (Rieke & Lebofski 1985), a value 
which agrees with that used in Babusiaux & Gilmore (2005) and 
also with that obtained by Indebetouw et al. (2005) by combin- 
ing data from the Spitzer Space Telescope and 2MASS. 

With the reddening-independent magnitudes, any compact 
structural feature in the CMD would appear as a "clump" of 
stars, as they will be located at approximately the same distance 
from the Sun and hence they would have the same apparent K- 
magnitude. As the red-clump stars are the dominant giant pop- 
ulation (Cohen et al. 2000; Hammersley et al. 2000), we can 
derive the distance to the feature by a non-linear least-squares 
fit of the sum of a Gaussian function and a second-order poly- 
nomial (eq. |2]i to the histogram of red giant stars (Stanek et al. 
1994; Stanek & Garnavich 1998; Babusiaux & Gilmore 2005): 



N(m) — a + bm + cm 2 ■ 



N RC 

CT RC V27T 



exp 



(ra RC - m) 1 



2 °lc 



(2) 



The Gaussian term corresponds to a fit of the red-clump pop- 
ulation of the bulge/bar, while the second-order polynomial re- 
flects the contribution of the background stars. The combined ef- 
fect of distance plus extinction produces an increase in the mix- 
ture between faint dwarfs and giants. In order to subtract the 
contribution of foreground dwarf stars we have used the SKY 
extinction model (Wainscoat et al. 1992) to derive the theoreti- 
cal traces in the diagrams for the giant population in each field. 
With these theoretical traces we can isolate the population of red- 
clump stars in each CMD; thus, the fit is made only after remov- 
ing the contribution of dwarf stars in the reddening-corrected 
counts (Fig.[T|). 

To transform the derived extinction-corrected magnitudes 
into distances from the Sun we use: 



Hk = K e + 



-(J - K s ) - M K . 



(3) 



To estimate an upper limit to the reddening-independent 
magnitude that can be analysed reliably we have followed the 
same procedure described in Babussiaux & Gilmore (2005). 
First, an estimate of the completeness limits in J and K (Jq, Kq) 
are derived by determining when the number counts as func- 
tion of magnitude stops increasing. The relative completeness of 
K e is then derived by combining these estimates with the mean 
colour of the giants in the field ( J — K)m, and adding a 0.5 mag 
safety margin to take into account the spread in the mean colour: 

max(K e ) = min(J c -(J- K) M , K c ) - [^ L < J ~ *0m) + 0.5.(4) 
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Fig. 1. Two examples of the use of the theoretical traces for the 
giant population by means of the SKY model. Stars isolated be- 
tween both traces (corresponding to the K0 and the M0 popu- 
lations) are assumed to be red-clump stars and are extracted to 
analyse their spatial distribution. 



Figure [2] shows all of the analysed fields used in this work. 
There are three different kind of fields: 

- Those (labelled "discarded" fields) in which the complete- 
ness made it impossible to fit the red-clump distribution ac- 
curately. 

- Those (labelled "rejected" fields) in which no successful fit 
was possible, as no clear structural component was present. 
In these cases, the counts increase continuously with the dis- 
tance modulus until reaching the completeness limit for the 
field. These in-plane fields are purely disc fields and are lo- 
cated too far above or below the galactic plane for either 
the bulge or the bar components to be observed, so that no 
"hump" in the red-clump counts is expected (Fig. [3j. 

- Those marked as "used" in Fig.|2]where it was possible to de- 
rive a distance to the red-clump population. After discarding 
all those fields which were not useful, we obtained results in 
49 of the total sample of 205 fields. 
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Summary of the 205 TCS-CAIN fields analysed in this 



1=27° b=2° 



200- 




12 13 14 
distance modulus 

Fig. 3. Example of a field where no peak in the counts of the 
red-clump distribution is observed, as we are observing a field 
where only the disc component is present, hence the star counts 
increase up to the completeness limit of the field (shown as a 
vertical dashed line). 



b < 5° were compared where JHK S photometry was available; 
also, there are no effects of completeness in either H e or K e . In 
general, the fits using the two different reddening-independent 
magnitudes are all consistent to within 0.2 mag, with a mean 
difference of 0.05 mag (cr = 0.11 mag) for a total of 18 fields 
(see Figs.|4]and|5]l. Hence, the derived distances using (J-H, H) 
data can be compared with those obtained using ( J - K, K) data. 



E 




Fig. 4. Comparison between the distance modulus obtained with 
//-band data with that obtained with /T-band data for the same 
fields. The solid line shows the 1:1 ratio. The two quantities are 
well correlated and give similar distances. 



4. Model fits 

The distances to the bar/bulge red-clump cluster of stars are sum- 
marized in Table Q] and examples of the fits are shown in Figure 
[6] Error estimates for the distances have been obtained as a com- 
bination of the systematic uncertainties in the intrinsic colour 
and absolute magnitude of the red-clump population and the fit- 
ting error of eq. [2] added in quadrature. The dispersion in dis- 
tance is represented in each case by cr. This is estimated from a 
simple deconvolution: 



3.1. Distances derived with H/J - H data 

Some fields were observed only in the J and H bands, since 
far from the Galactic plane the counts in the K and H bands 
are nearly equivalent due to the low extinction. For this reason 
some of the fields in this work use //-band extinction-corrected 
magnitudes (H e ): 

H e - H - Ah (J-H), (5) 

fi H =H e + A " ( J - H) - M H , (6) 
Aj-A H 

using (J - H) = 0.5, Mh = —1.5 (Wainscoat et al. 1992) and 
A H IE j-h = 1.41, following Rieke & Lebofsky (1985) and in 
agreement with the result of Nishiyama et al. (2006a). These 
values were used in eight fields, those with \b\ = 3.5° and \b\ = 
7.5°, and three additional fields at b = 5° (those at I < 0°). 

To test the reliability of deriving distances using either K- or 
//-band corrected magnitudes, the results for the fields with 2° < 



where cr RC is the dispersion of the Gaussian fitted to the his- 
tograms by eq. [2] cr is the intrinsic dispersion of the red-clump 
luminosity, which has been estimated in the range 0.15-0.2 mag 
(Alves 2000; Lopez-Corredoira et al. 2002) and finally, cr e is the 
photometric error, which for TCS-CAIN is of the order ~0.08 
mag (Cabrera-Lavers et al. 2006). 

4.1. Contamination from the dwarfs 

Although we have used the theoretical traces from the SKY 
model to isolate the red-clump population to avoid contamina- 
tion due to other populations in the counts (Fig. [TJ, there will 
be some contamination from other sources and the effect of this 
on the overall shape of the reddening-corrected counts should be 
estimated. For this, we have used the SKY model to predict the 
relative number of dwarf stars that will fall into the red-clump 
region of the CMDs (between the dashed lines in Figure[T]i- The 
dwarfs predicted by the SKY model that fell into the red-clump 
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Fig. 5. Spatial distribution of red-clump stars obtained with ei- 
ther H- or ^T-band corrected magnitudes for the fields I = 7°, 
b = 5° (above) and I = 5°, b = 5° (below). The distributions 
are very similar, and the maxima are nearly coincident. Vertical 
lines show the limiting magnitude derived using the K (solid) 
and H (dashed) bands, respectively. 
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Fig. 6. Histograms of the distance modulus in three of the fields 
used in this work. Fits of eq. [2] to the histograms are shown as 
solid curves, while vertical lines show the limiting magnitude 
estimates up to which the fits were performed. 



area on the CMD were subtracted from the measured counts to 
give a new "corrected" distribution. A full description of this 
procedure can be found in section 3.3.3 of Lopez-Corredoira et 
al. (2002). What is found is that for m K < 12.5 only 2.5-5% of 
the detected sources are dwarfs, but that this rises to 10—40% for 
13 < m K < 14. 

As shown in Fig. |7j both for an in-plane and an off-plane 
field, there is no significant difference in between the corrected 
and uncorrected distribution. The effect of dwarfs is more im- 
portant for fainter magnitudes (hence for apparently more distant 
sources), so the distributions overestimate the number of sources 
for higher values of the distance modulus. However, the peak lo- 
cation is dominated by the red-clump population, and any possi- 
ble contamination does not affect its location. We fitted eq.|2]to 
the "corrected" distributions and obtained a mean difference in 
the distance modulus derived from either the uncorrected or cor- 
rected data of 0.05 mag (which, translated into distances, yields 
to a difference of around 100 pc in this case), which is lower than 
the accuracy of the fit itself. This calculation was done only for 
a test sample of four randomly selected fields. Since we did not 
find any significant difference in the results, correction was not 
attempted for with the remaining fields. 



Lopez-Corredoira et al. (2002) demonstrated that the effect 
of contamination is more severe off the plane, as extinction helps 
to separate the populations. The results for off-plane fields show 
a much better correction of the counts, giving a slightly narrower 
distribution around the peak. Although applying this could in 
theory improve the fit, the effect is small and introduces another 
factor based on models, hence we decided not to apply this cor- 
rection. 

5. The Distribution of sources in the inner Galaxy 

5.1. In-plane fields (\b\ < 1°) 

Results for in-plane fields in the range 18° < I < 29° show a 
narrow distribution of distances from the Sun. This is consis- 
tent with a high-density feature seen in the red-clump stars with 
position angle 40?5 ± 3°.9 (after a linear fit through the points) 
with respect to the Sun-Galactic Centre direction. First panel in 
Figure[8]shows the distribution of red-clump density maxima for 
a face-on view of the Galaxy, with the Galactic Centre at (0,8000 
pc). The solid line shows a position angle of 45° (consistent with 
the 43 degrees given for the bar in Hammersley et al. 2000). The 
dashed line shows a position angle of 25°, which is consistent 
with the results of Dwek et al. (1995), Lopez-Corredoira et al. 
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Table 1. Distance modulii and distances derived with the red- 
clump stars for the 49 fields used along this work, with each 
field covering 0.07 deg 2 on the sky centred on the Galactic coor- 
dinates given in Cols. (1) and (2). The last eight fields are those 
in which //-band data were used to derive distances. 



1=28 b=o 



m 


b(°) 


fi (mag) 


D (kpc) 


<x (kpc) 


18 


0.0 


14.01±0.07 


6.317+0.229 


0.47 


20 


0.0 


13.68±0.07 


5.449+0.199 


0.44 


21 


0.0 


14.05±0.09 


6.455+0.298 


0.51 


22 


0.0 


13.95±0.08 


6.187+0.238 


0.47 


26 


0.0 


13.81+0.07 


5.778+0.186 


0.41 


27 


0.0 


13.71+0.08 


5.525+0.216 


0.42 


28 


0.0 


13.66+0.11 


5.392+0.303 


0.52 


29 


0.0 


13.62+0.10 


5.308+0.258 


0.50 


20 


-0.5 


13.36+0.08 


4.704+0.184 


0.65 


20 


0.5 


13.79+0.08 


5.719+0.221 


0.51 


26 


-0.5 


13.58+0.08 


5.211+0.204 


0.52 


26 


0.5 


13.65+0.12 


5.366+0.304 


0.54 


27 


-0.5 


13.55+0.08 


5.126+0.202 


0.56 


27 


0.5 
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0.43 


7 


5.0 


14.04+0.09 


6.430+0.271 


0.72 


-5 


5.0 


14.67+0.11 


8.588+0.459 


0.43 


-3 


5.0 


14.52+0.09 


8.015+0.340 


0.49 


-2 


5.0 


14.56+0.14 
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0.44 
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(2005), Stanek et al. (1997) and, more recently, Babusiaux & 
Gilmore (2005). The range of longitudes covered in these fields 
(18° < I < 29°), however, is very different from those which 
obtain the smaller position angle. Only a single point at I = 20°, 
b = -0.5°, is observed near the latter position angle. This point 
could be due to local highly anomalous extinction (see Sect. 15. 31 ), 
as this line of sight does run close to a major star formation re- 
gion (Garay et al. 1998; Sridharan et al. 2005). Furthermore, the 
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Fig. 7. Histograms of the distance modulus for two of the ob- 
served fields. Points/solid lines show the histograms derived di- 
rectly from the data, while dotted/dashed lines show the distri- 
butions once corrected for dwarf contamination (using the pre- 
dictions of the SKY model). Vertical dashed lines stand for the 
completeness-induced limits for each field. 



over-density is is apparently somewhat inhomogeneous, as was 
noted by Picaud et al. (2003) at I = 20° and I = 21°. 

The gap in the longitude range 23° < I < 25° is due to the 
very high extinction that extends up to \b\ ~ 0.4°, which is no- 
ticeable also in the mid-infrared range (see Fig. 1 of Benjamin et 
al. 2005). This makes it impossible to obtain an accurate fit for 
the fields located in this region. The same problem applies to the 
in-plane data for 0° < I < 15°. In these fields the completeness 
limit is far brighter than the apparent magnitude of the K gi- 
ants for any structure that would be present, that is significantly 
fainter than the completeness limit. Hence, the above method 
cannot be applied. However, in the CMD the giant branch from 
the feature is still visible but its distance cannot be measured. 



5.2. Off-plane fields (\b\ > 1°) 

The results for \b\ = 1° (second panel in Fig. [8]) show that there 
are two distance regimes for the red-clump population, one that 
follows approximately the 43°4 + 3°5 angle (nine points) and a 
second with a smaller angle of 21°9 ± 1°0 (six points). Most 
of the points in this latter feature are at / < 15°, while the 
points associated with the former are mostly located at I > 20°. 
Nishiyama et al. (2005) also observed the red-clump stars at 
\b\ = 1° over the range 8° < I < 15° and there is good agree- 
ment between their distribution of red-clump maxima and ours 
(Fig. [9]l, even in the central regions of the Galaxy. 
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Fig. 8. Spatial distribution of red-clump giants maxima in the 
.XT-plane at different latitude intervals, with the Sun at (0,0) 
and the Galactic Center at (0,8 kpc) marked with a big circle. 
Two possible configurations for the observed distribution are 
also shown: a feature with a position angle of 45° (solid line), 
and another one with a position angle of 25° (dashed line). Dot- 
dashed lines define two circles with radii 4.5 kpc and 3 kpc, re- 
spectively, while different lines of sight towards the inner Galaxy 
are also shown in intervals of 10° in longitude. Error bars have 
been estimated from the distance uncertainties, assuming these 
divided equally in both the X- and F-axes. 
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Fig. 9. Spatial distribution of red-clump giants maxima pre- 
sented face on for data at \b\ = 1° and the results obtained for 
Nishiyama et al. (2005) for comparison (shown as open circles). 



For \b\ - 2° the results are less conclusive (third panel in Fig. 
[§). This is because we are not deep enough to be complete in 
the red-clump giants in the innermost Galaxy (|/| < 5°) and this 
makes it difficult to observe the peak in the red-clump counts. 
Also, for fields \l\ > 15° there is no red-clump cluster observed 
in the modulus distance histograms, i.e. there is no feature there. 
Only one point is detected with I > 15° and this has the highest 
dispersion in distance of any field. The remainder of the fields lie 
where there is little difference in distance between the two po- 
sition angles, making it difficult to clearly assign the red-clump 
density maxima to one of them. A linear fit to those data yields 
a global position angle of 36°6 ± 5°8 (after discarding those at 
/ = 17°), although if the points are assigned to the position angle 
that runs closer to that position, and are then averaged separately 
we obtain position angles of 44°. 1 + 2°8 and 25°. 1 + 2°7, respec- 
tively. 

Farther away from the plane {\b\ > 3°) it is possible to 
sample the innermost Galaxy as the incompleteness effect dis- 
appears; however, there are no observable features in the his- 
tograms for I > 1° as the disc component dominates and there 
is no dense structure at a specific distance. The distribution of 
red-clump peaks shows in this case a position angle of 30°7 + 
5°2 with respect to the Sun-Galactic Centre line (lower panel 
in Fig. [8j. These results are compatible with the position angle 
of 29° ± 8° derived in Lopez-Corredoira et al. (2005) for the 
bulge of the Milky Way using 2MASS data. Also, this result is 
consistent with the position angle of 22°+5°5 recently derived 
by Babusiaux & Gilmore analysing their private in-plane data at 
I/I < 10°. 



5.3. The effect of extinction on the results 

A global extinction law is assumed when deriving the reddening- 
independent magnitudes. However, recent results show that there 
are noticeable differences in the extinction law derived in differ- 
ent directions towards the inner Galaxy (Nishiyama et al. 2006a) 
and even very small variations in the ratio of total to selective ex- 
tinction are present on very small scales (Gosling et al. 2006), so 
care must be taken when applying reddening corrections. The 
combination of this plus the possible mixture of sources com- 
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Table 2. Fitted position angles and mean dispersion in distances 
for the different latitude ranges. 



b range 


N 


0bar (°) 


0bulge (°) 


(o-) (kpc) 


\b\ < 1° 


14 


40.5±3.9 




0.48+0.04 


\b\ = 1° 


15 


43.4+3.5 


21.9+1.0 


0.55+0.09 


\b\ = 2" 


7 


44.1±2.8 


25.1+2.7 


0.51+0.03 


\b\ > 3.5° 


13 




30.7+5.2 


0.45+0.02 


|*| = 3.5° 


3 




21.6+3.9 




\b\ = 5° 


8 




31.5+3.5 





ing from the two position angles would broaden in the distance 
histograms, hence giving a higher dispersion in distances, cr. 

To examine this effect, we have calculated the mean disper- 
sion in distances in the four latitude ranges showed in Fig. [8] 
(Table |2j. There is clearly less dispersion at \b\ < 1° and at 
\b\ > 3° than at intermediate latitudes, where the mixture of 
sources becomes more important (error bars shown in Table [2] 
account only for the standard deviation of values). The mean 
values all cluster around 0.5 kpc, but again, with slightly larger 
values at \b\ = 1° and \b\ = 2°. Hence, if we consider the re- 
sults only for those fields that present <x < 0.5 kpc we obtain the 
red-clump peak distributions shown in Figure [TOl 

It is noticeable that in all the cases those points that were 
originally located at positions intermediate between the two po- 
sition angles have been removed. Furthermore, only one point 
has disappeared from the fields at \b\ < 1° and \b\ > 3° as a 
reflection of the lower dispersion values obtained in those lati- 
tude ranges. The point removed is that noted in ^5. H at I = 20°, 
b = -0.5° suggesting that extinction is responsible for its dis- 
crepant location. 

5.4. Discussion of the two position angles 
The above shows that: 

- The red-clump stars show that there is a dense feature seen 
on the plane in lines of sight 0° < I < 28°, but at \b\ > 2° it 
is only seen for 0° < / < 10°. In other locations there is no 
dense feature, and only the disc is seen. 

- Away from the Galactic plane for \l\ < 10° the peaks in den- 
sity of the red-clump population follow a very well defined 
feature with a position angle of 23°. 1 ± 0°9 with respect to 
the Sun-Galactic Centre line. 

- In the plane (\b\ < 1°) for \l\ > 18° there is a longer structure 
with a distinct position angle of 43°0 ± 1°9. This feature is 
not seen more than 2 degrees away from the plane. 

- At intermediate latitudes, there is a mixture of sources com- 
ing from both structures producing a more spread out distri- 
bution in the distances derived from the red-clump location 
in the CMDs. 

The position angles are derived by a weighted mean of the values 
summarized in Table|2]for both features. 

Sevenster et al. (1999) found a similar effect when compar- 
ing A^-body models with their sample of OH/IR stars and noted 
that the position angle changed from <p = 44° when using in- 
plane fields to <p — 25° when only high latitude fields were con- 
sidered. They also noted in the literature that lower values for the 
viewing angles determined from stellar data were found when 
low-latitude fields were excluded from the fit, which is in agree- 
ment with the work presented here. 
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Fig. 10. Same as in Fig. [8] but considering only those fields 
where dispersion in distances is less than 0.5 kpc. 



The feature seen on the plane 18° < I < 21° in the red-clump 
stars can have no direct link with the Scutum spiral arm. The 
tangential point to the stars in this arm is at I ~ 32°, and it would 
be expected that any features associated with an arm would be 
strongest there. However, there is no red-clump cluster seen at 
that location (Hammersley et al. 2000), and inwards of the this 
point the effect of the arm would become less (Hammersley et al. 



A. Cabrera-Lavers et al.: Tracing the long bar with red-clump giants 



9 



1994). Furthermore, a spiral arm would not form a feature that 
appeared to run directly towards the Galactic Centre. 

Extinction can also be discounted. There are a range of ex- 
tinctions covered between 18° < I < 29° and in no field is the 
extinction less than is predicted by the simple extinction model 
used in Wainscoat et al. (1992). In any case, the method used 
corrects for extinction. Furthermore, our result is in extremely 
good agreement with the result of Benjamin et al. (2005) (see 
Fig.fTTTi. where the effects of extinction are significantly reduced. 
Also, in Fig.QT]we show how the concentration of data seen at 
18° < I < 27° is not an isolated feature, and how it is consistent 
with other estimates for this structure obtained with deeper data 
than ours in different lines of sight to the inner Galaxy. In fact, 
it is clear that this feature runs to the inner Galaxy, at least up 
to 1.5-2 kpc from the Galactic Centre, where it merges with the 
bulge. 




-2000 
X(pc) 

Fig. 11. Spatial distribution of red-clump giants maxima for in- 
plane fields in the .XT-plane, with some of the more recent esti- 
mates for this feature (H00: Hammersley et al. 2000; P03: Picaud 
et al. 2003; B05: Benjamin et al. 2005). Note there is an excel- 
lent agreement between the different sets of data. The long bar 
with a position angle of 43°0 derived in this work is represented 
by a solid line (and a width of 1 kpc is represented by dashed 
lines). 



The only scenario that readily fits the star count data and 
CMDs, those presented both here and elsewhere, is that of a 
triaxial bulge plus a long bar. The bulge has a position angle 
of about 25° and dominates the counts to |/| ~ 12° and up to 
\b\ = 7° or farther off the plane. The in-plane bar has a posi- 
tion angle of 43° with a half-length of some 4 kpc but is only 
detected close to the plane. This bulge + bar picture is not so 
unreasonable. Recent A^-body simulations on the secular evolu- 
tion of disc galaxies predict that in a bar-driven evolution, the 
bar has a thick inner part of shorter extent and a thin outer part 
of larger extent (Athanassoula 2006), a fact that is in good agree- 
ment with observations (Athanassoula 2005; Bureau et al. 2006). 
Also, the ratio of the thin bar size to the inner bulge length that is 
derived from our data (between 2.2-2.8 assuming a bulge length 
of 1.5-2 kpc) is in completely agreement with that was observed 
in external galaxies (2.7 + 0.3, Liitticke et al. 2000). 

In order to illustrate this, Fig. [12] shows the results from 
Babusiaux & Gilmore (2005; hereafter BG05). Also shown on 
the plot are the results from this work and position angles of 



45° and 25° degrees. The BG05 data set is deeper than ours (as 
a result of a pixel scale five times smaller), permitting them to 
reach the inner Galaxy with in-plane data. It is noticeable that the 

BG05 points at / ~ 10° and I 5° coincide with positions for 

the triaxial bulge obtained here; however, the point at I 10° 

lies on the 45° position angle and is well away from the position 
required for the triaxial bulge. 




-4000 -2000 
X(pc) 



2000 



4000 



Fig. 12. Distribution of positions for the red-clump sources ob- 
tained in this work at \b\ < 1° and \b\ > 3° compared with the 
results of Babusiaux & Gilmore (2005) by using in-plane NIR 
data. Dot dashed lines delimite circles with radii 4.3 kpc, 3 kpc 
and 2 kpc, respectively. 



5.5. Scale height of bar/bulge sources 

From the histograms of the red-clump sources an estimate can 
be made of the mean scale height of the sources associated with 
the long thin bar. We have used only those fields with / > 20° 
and \b\ < 1°, as at b — 2° there are not sufficient bar sources 
for a reliable extraction. For each field the space densities along 
different lines of sight were determined by selecting the stars 
within one cr of the maximum of the histograms of the distance 
modulus (see Table[T). Combining the distances along the line of 
sight to the features (d) and the Galactic coordinates of the fields 
(I, b) allows these space densities to be expressed with respect to 
the galactocentric distance (R) and the height (z) above the plane, 
taking into account the location of the Sun 15 pc above the mean 
plane (Hammersley et al. 1995). The result is shown in Figure 
fT3l It can be observed that the mid-point of the distribution is at 
Z = -41 pc. Hence, we have represented in Figure [14] the space 
density with respect to a vertical distance to the plane z*, defined 
as z* = \z + 41 1. This means that the location of this hypothetical 
long thin bar is ~ 0.5° tilted with respect to the mean plane, 
a result that is compatible with the values found in Picaud & 
Robin (2004) for the orientation of the outer bulge. 

We have also tested a sech 2 profile to compare with the ob- 
served space density for a simple exponential. The two equations 
used are: 

p = po exp(-z*//i z ), 
p = posech 2 (z* /h z ) 

where h z is the scale height of the bar sources and po is the space 
density in the Galactic plane. The exponential gives h z = 99.3 + 
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Fig. 13. Variation of the space density as a function of distance 
from the Galactic plane. For comparison, the best fit of an expo- 
nential law (solid line) and a sech 2 law (dashed line) are shown. 
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Fig. 14. Vertical dependence of the density of bar stars. The best 
exponential fit, with a scale height of 99.28 + 3.13 pc, is also 
shown as a solid line. 

3.1 pc and po = 0.022 + 0.014 pc~ 3 , whereas the sech 2 profile 
yields h z = 105.4 + 5.3 pc and p = 0.018 + 0.013 pc~ 3 . Both fits 
are shown in Figure [T3l 

An exponential scale height for the bar of around 100 pc 
is a factor 2 higher than that previously obtained by Lopez- 
Corredoira et al. (2001) or Hammersley et al. (1995) from NIR 
star counts. However, it should be noted that here we are look- 
ing at a very different star type than in the previous studies. In 
both of those cases the stars observed were considerably more 
luminous, and spectra of a selection of these sources (Garzon 
et al. 1997) showed that there were a large number of young 
sources at I = 27° and I = 21°. The red-clump sources observed 
here probably represent an older population and hence the larger 
scale height would be expected. However, the scale height for 
the red-clump stars in other Galactic components is far higher. 
For the thin disc it is ~260 pc (Cabrera-Lavers et al. 2005) and 
applying the same method used above to the bulge for the fields 
at \b\ > 3.5° gives a scale height of h z = 488 ± 28 pc (see Fig. 
[151 ). Hence the scale height is a further clear indication that on 
the plane between 18° < / < 29° there is a geometric component 
that is very different to the bulge or thin disc. 




|z| (pc) 



Fig. 15. Vertical dependence of the density of inner bulge stars. 
The best exponential fit shown as a solid line now has a scale 
height ~5 times higher than that of the thin bar. 

6. Systematic uncertainties in the red-clump 
method 

6. 1 . Deriving distances from the number count histograms 

The distances to the red-clump population were obtained by fit- 
ting eq. [2] to the distance modulus histograms. However, the 
maximum of star counts vs. magnitude is not strictly coincident 
with the maximum in density along the line of sight, which is 
proportional to the first multiplied by a factor 1/d 3 . Hence the 
magnitude histogram counts, N(m), should be transformed into 
density along the line of sight, p(d), and eq. [2] then fitted to the 
resultant distribution. From the relationship between p(d) and 
N(m) it can be shown that the difference between the corrected 
distance to the maximum in the density distribution (r* n ) and the 
distance obtained to the maximum in the counts histograms (r m ) 
is: 

r* m = r„, + Ar m (8) 

= w 

r m p (r m ) 

r m being the distance to the maximum in the N(m) histogram 
(hence, N'(m[r m ]) = 0) and p"(r m ) the second derivative of p 
with respect to the distance, measured at r m . As p"(r m ) < 0, 
then r* m < r m , thus the corrected distances are slightly lower than 
those obtained from the maximum of the magnitude histograms. 

This effect is more noticeable for the bulge component than 
for the long thin bar, as p"(r m ) is very large for this latter. We 
have estimated the range of values for Ar m for the bulge and 
for the long bar, obtaining the true density distributions along 
the line of sight. We have used those fields with \b\ > 3.5° as 
representative of the bulge, while those in-plane fields with / > 
20° are representative of the long bar. It is found that Ar,„ is in the 
range 100-300 pc in bulge fields, whereas for long bar fields the 
effect is far smaller, in the order of 25-50 pc. Figure [161 shows 
two examples of the change in the overall shape of the density 
distribution compared to the number count histograms used to 
derive the distances to the red-clump population. The effect on 
the determined geometries for either the bulge and the long bar 
is not large in either case. For the long bar it is ten times lower 
than the uncertainties in the fit itself (see Table [TJ and even for 
the bulge it is only about 50% of the uncertainty in the fit. 
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6.2. Difference in maximum density along the line of sight 
and the position of the major axis of a triaxial structure 

The position of the maximum density along the line of sight is 
not coincident with the the position of the major axis of a triaxial 
structure unless the structure is very thin along the line of sight 
(a bar). There are two effects: 

- The more important one is that the density along the line 
of sight reaches a maximum at the tangential point of the 
innermost ellipsoid along it and this is not in general on the 
major axis (see Fig. \TT}. 

- For off-plane regions, the lines of sight are not parallel to the 
b — 0° and so the larger the distance the greater the height 
above the plane, and hence the lower the density. Therefore, 
the maximum density along the line of sight is closer than 
the real maximum in a plane parallel to b — 0° . 
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Fig. 17. Schematic representation of the difference between the 
maximum density position along the line of sight (filled square) 
for triaxial ellipsoids and the intersection of their major axis with 
the line of sight (filled circle). 



that in the n=2 case corresponds to isodensity ellipsoids. A and 
B are the axial ratios of the second and the third axes with re- 
spect to the major axis of the ellipsoids, and x, are the Cartesian 
coordinates with the axis of the ellipsoids concentric with the 
centre of the Galaxy (see Fig. [TTT l. We assume that the minor 
axis is perpendicular to the Galactic plane (x$ = z). x, y, z are the 
Cartesian coordinates with XY defining the plane of the Galaxy 
with the y-axis in the Sun-Galactic Centre line and a is the angle 
between the major axis of the ellipsoid and this y-axis (Fig. [17}. 

The tangential point of the line of sight with the maximum 
density (filled square in Fig. [TTTl. follows: 



dt 

T-(r m ) = 0. 
or 

Equation[TT] together with Eq. [TO] lead to 



A l (r m )^-(r m ) + 



4~ 1 (.r m )dx 2 



j"'(rm) dx 3 
B" dr 



8r A" 

which finally lead to 



dr 



(r m ) + - 



(ID 



(r m ) = 0,(12) 



n-l 



Fig. 16. Comparison between the density (solid line) and the 
number count (dashed line) distributions with respect to the dis- a > r 'm = 

tance along the line of sight for long bar (above) and bulge (be- 
low) fields as examples. Density distributions are arbitrarily dis 
placed along the F-axis in order to make both distributions com 
parable. 



(13) 



/■=o 
where 



We can derive analytically the difference (Ar = r m - r a ) 
between the maximum density position along the line of sight 
(r m ) and the intersection of the line of sight with the major axis 
(r a ) for a triaxial structure with monotonically density decreas- 
ing outwards. The isodensity contours are defined by the points 
in space with the same value of f: 



tan(/ + a) sin" X {1 + a) 
A" 



a„_j =(-1)" cos"~ 1 (/ + or) + 



tan" b 



B" cos{l + a) 



(14) 



and 



(xl + (x 2 /AT + (x 3 /BT) 



l/n 
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X^cos" 1 ' ffcos'(7 + a) + 
V/e [0,71-2], 



sin" 1 ' a sin' (7 + a) tan(7 + a) 
A" 



(15) 



while the position of the major axis (simple geometry with 
the application of the sine rule, see Fig.[T7li is 



R 



cos b sin(Z + a) 
In the simplest case, for n — 2, Eq.[T3l yields 



coso L 



cos a + 



tan(/+a-) sin a 
A 2 



cos(/ + a) + jt sin(l + a) tm(l + a) + (^) 



l 



(16) 



(17) 



COS(/+Q') 



Both of the expressions for r m and r a are coincident for A <K 
(very elongated ellipsoids) and tan b <sc (in the plane), but other 
cases are affected by a significant systematic error Ar. Two tri- 
axial structures were examined in order to estimate this effect: 

- A - 0.5, B = 0.4, a = 28°, typical of the bulge (Lopez- 
Corredoira et al. 2005): Fig. HI 

- A = 0.11, B = 0.04, a = 43°, typical of a long bar: Fig.[19j 
The low depth in the line of sight of the bar (around 1 kpc) 
is justifiable because of the low dispersion of the red-clump 
giants. 

The result is that Ar is negligible for the bar, with less than 
~150 pc of systematic error; however, it is not negligible for the 
bulge, which reaches a discrepancy of up to 1000 pc. 




Kdes.l 



Fig. 19. Difference between the maximum density position along 
the line of sight for triaxial ellipsoids of axial ratios 1:0.1 1:0.04, 
orientated with respect to the Sun-Galactic Centre line at an an- 
gle of 43° (expected for the long bar) and the intersection of their 
major axes with the line of sight. 



a = 13°3^ 
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The result for a is notably smaller than that directly de- 
rived from the fit to the red-clump density peak, but it is consis- 
tent with other estimates found in previous works (Freudenreich 
1998; Lepine & Leroy 2000; Picaud & Robin 2004). The axial 
ratios of the bulge are more or less coincident with those con- 
sidered in Lopez-Corredoira et al. (2005) and others (e.g. Stanek 
et al. 1997; Bissantz & Gerhard 2002; Picaud & Robin 2004). It 
should be noted, however, that eq. [17] is valid only for a series 
of triaxial ellipsoids, and so will not be correct if the bulge has a 
different geometry. However, when we use the general formula 
given by eq. [13] the results do not change noticeably (see Table 
0, all of them being compatible with a structure with a position 
angle of 10-15° (with a weighted mean of 12°6 + 3°2). 

Table 3. Parameters derived for different assumed bulge geome- 
tries 



Fig. 18. Difference between the maximum density position along 
the line of sight for triaxial ellipsoids of axial ratios 1:0.5:0.4, 
oriented with respect to the Sun-Galactic Centre at an angle of 
28° (expected for the thick bulge) and the intersection of their 
major axes with the line of sight. 
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A 


B 


«C) 




Type 


2 
3 
4 


0.48-- 
0.50±0.05 
0.55±0.03 


0.28^ 
0.45±0.05 
0.45±0.03 
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To avoid this effect in deriving the real parameters of the 
bulge we performed a chi-squared fit to the observed distances of 
the red-clump density maxima by means of eq. Q7)> which con- 
tains only three free parameters (A, B and a). We have selected 
only those points at \b\ > 3°, as they are more representative of 
this component. For each set of model parameters we calculate 
X 1 - Once this was calculated for all sets in the parameter space, 
we select that one which gives a minimum^ 2 . Again, for n — 2 
the values which give the minimum^ 2 are: 



In summary, the position angle we have derived in this work 
for the bulge (73°. 1 ± 0°9) by using the distribution of red-clump 
density maxima might be overestimating the real position angle 
by around 8-10°. So the real position angle of the inner bulge 
is around 13-15° with respect to the Sun-Galactic Centre line 
(nearly coincident with the value of 12°6 ± 3°2 obtained via 
chi-squared fitting). In the case of the long thin bar, the effect 
analysed in this section is less significant and so do not alter the 
geometry of this component. 
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7. Confirming the result using 2MASS star counts 

The 2MASS survey (Skrustkie et al. 2006) provides a complete 
coverage of the Galactic Plane in the JHK S . Unfortunately, the 
2MASS data are not deep enough to reach the red-clump stars in 
the inner Galaxy, hence the 2MASS data cannot be analysed us- 
ing the above method. However, the 2MASS star counts can be 
useful to confirm the proposed morphologies in the inner Galaxy, 
as they are a powerful tool when searching for asymmetries in 
the stellar distribution. While the bulge is observed in off-plane 
regions up to b ~ 10° and at |/| < 15°, the long bar is visible 
only in the in-plane regions, \b\ < 2° and up to / = 27° at posi- 
tive longitudes and but only to / = -15° at negative longitudes 
(from the length and position angle). For this study, counts up 
to trig = 9 mag will be used as this gives a high inner Galaxy- 
to-disc contrast (Garzon et al. 1993) whilst still providing suffi- 
cient sources to give good statistics. A similar analysis was per- 
formed in Lopez-Corredoira et al. (2001) with the combination 
of DENIS (Epchtein et al. 1997) and TMGS (Garzon et al. 1993) 
/f-band counts. The advantage of using 2MASS is that it covers 
the whole region homogeneously and so makes testing the model 
simpler. 

Figure [20] shows the 2MASS star counts with nix < 9 
mag (averaged over Al = 5°) at three different latitude slices: 
-0.25° < b < 0.25°, -1.25° < b < -0.75° and -2.25° < 
b < -1.75°. At the distance of the Galactic Centre (7.9 kpc), 
the heights above the Galactic plane are pc, -140 pc and -280 
pc respectively. 

At b — 0° there is a large-scale asymmetry about / = 0° 
with far more counts at positive longitudes than at negative lon- 
gitudes. Farther from the plane young inner Galaxy features with 
very small scale heights should not be significant and the ef- 
fect of extinction is reduced. Hence, although the contrast with 
the disc is reduced there should still be sufficient bar counts to 
see something and the accuracy of model should be massively 
improved as it cannot deal with the patchy extinction on the 
plane. In Lopez-Corredoira et al. (2001) it was shown that inner 
Galactic star counts are more or less symmetric in latitude, with 
small differences attributed to the Sun being about 15 pc above 
the Galactic plane (Hammersley et al. 1995) and differences in 
the extinction above and below the plane. Here, we concentrate 
our analysis only on those slices at negative latitudes and we 
assume that the structure of the inner Galaxy is more or less 
symmetric about the Galactic plane. 

We have constructed a very simple Galactic model that in- 
cludes the contributions of the Galactic disc, a bulge, and a long 
thin bar to reproduce the 2MASS counts at {b) = -2.0°. For the 
Galactic disc we have used the luminosity function of Eaton et 
al. (1984) and an exponential density for the outer Galaxy with 
a scale height of 285 pc and a scale length of 2.1 kpc (Lopez- 
Corredoira et al. 2002). The inner 4 kpc of the disc, however, has 
a constant density (Lopez-Corredoira et al. 2004). It has has been 
shown that continuing the exponential disc to the Galactic Centre 
leads to a significant overestimation of the TT-band star counts 
(Lopez-Corredoira et al. 2001; Benjamin et al. 2005). This result 
is consistent with the disc being a Freeman type II disc, which 
are common in galaxies with long bars. The triaxial bulge has 
ratios 1:0.5:0.4 and a position angle of the major axis with re- 
spect to the Sun-Galactic Centre line of 27° (Lopez-Corredoira 
et al. 2005). The bulge luminosity function and density law have 
been also taken from this work. Finally, a simple bar model in 
agreement with Hammersley et al. (2000) was added. The bar 
has a depth along the line of sight of 500 pc, a half-length of 4 
kpc and a position angle of 43°. The distribution was assumed 
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Fig. 20. 2MASS star counts with Mk < 9 mag in three different 
latitude ranges. Counts have been averaged over Al = 5° to ob- 
tain a large number stars in each region, in this way minimizing 
the Poissonian error. 



to be constant along the bar, but exponential in height above the 
plane and the luminosity function used was the same as for the 
disc although the density was then normalized to make the total 
counts match those at I = 27°. In all the cases, the extinction 
model used is as described in Wainscoat et al. (1992). 

Whilst the model is simple, it can be shown from Fig. 1211 that 
it does reproduce the counts fairly well at b — -2°. The model 
predicts: 

- The increasing counts with decreasing longitudes for |/| > 
30° due to the exponential disc 

- The sharp jump in the counts at / = 27° due to the bar 

- The basically flat region between 12° < |/| < 27° due to 
exponential disc not continuing into the centre and the bar 

- The steep rise inwards of \l\ < 10° as the bulge dominates. 
The bar/bulge ratio in this region is below 0.025 in this re- 
gion 
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There are, however, some other features not well predicted by 
the simple model. As no ring was included in the model there 
are excess counts for the region -30° < / < -22°. Sevenster et 
al. (1999) also observed an excess of OH/IR stars at I = -22°. 
The possible implications for an elliptical ring are described in 
section 6. 1 of Lopez-Corredoira et al. (2001) but are beyond the 
scope of this paper, so will not be repeated here. Similarly, the 
model does not include the Scutum spiral arm so near the tan- 
gential point (/ = 33°) there could also be problems. A fuller 
analysis of the inner Galactic star counts can be found in Lopez- 
Corredoira et al. (2001; section 3). Here, we are only interested 
in testing the plausibility of the bulge + thin bar scenario when 
reproducing the observed counts. 




500 



40 30 20 10 -10 -20 -30 -40 
galactic longitude (degrees) 

Fig. 21. 2MASS star counts to K < 9 mag and the predictions of 
a model that includes the contribution of disc (solid line), bulge 
(dot-dashed line) and a long thin bar (dashed line). The sum of 
all the components (strong solid line) reproduces the observed 
counts very accurately. 

It is, however, clear from the above that a disc + thin bar 
+ triaxial bulge can account for the a large part of the observed 
counts, in particular all of the large-scale Galactic features. Were 
the bar not to be included, the fit would not be as good; in par- 
ticular the asymmetry in longitude would not be reproduced. 

8. Conclusions 

By analysing the distribution of red-clump sources along dif- 
ferent lines of sight towards the inner Galaxy, we have demon- 
strated that the Milky Way has two different structures coexisting 
in the inner R < 5 kpc: a long thin bar with a half length of 4 
kpc constrained to the plane and a position angle of 43°, and a 
thicker bulge with a position angle of 13-15°. 

In the TCS-CAIN data, each structure dominates at different 
ranges of longitude and latitude. At \b\ > 3.5°, \l\ < 10° the bulge 
is the only component that is observed whereas at 20° < / < 28° 
and \b\ < 0.5° the long thin bar dominates. The scale heights 
of these structures are also very different. The red-clump giants 
of long -bar sources have a scale height of ~ 100 pc, whereas the 
disc is over 2.6 times this thick and the bulge sources have a scale 
height ~5 times larger still. Hence the long bar is geometrically 
a different structure to the triaxial bulge or disc, although this 
does not necessarily mean that they are not related. 

We have also checked how there are some systematic un- 
certainties in deriving distances from the red-clump population 



peaks that have not been taken into account in other similar 
analyses (e.g. Stanek et al. 1994, 1997; Nishiyama et al. 2005; 
Babusiaux & Gilmore 2005). Those uncertainties make the po- 
sition angles derived in those works around 10° larger than the 
real ones. 

We have shown that the large-scale 2MASS star counts in the 
inner Galaxy can be accurately described by a combination of a 
triaxial bulge, a bar and a disc. Hence, discussion as to the angle 
of the bar should be separated into two separate discussions: the 
form of the triaxial bulge and the form of the long bar. 
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